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Abstract. We show that the multivariate generating function of appropri- 
ately normalized moments of a measure with homogeneous polynomial density 
supported on a compact polytope V C M'^ is a rational function. Its denom- 
inator is the product of linear forms dual to the vertices of V raised to the 
power equal to the degree of the density function. Using this, we solve the in- 
verse moment problem for the set of, not necessarily convex, polytopes having 
a given set S of vertices. Under a weak non-degeneracy assumption we also 
show that the uniform measure supported on any such polytope is a linear 
combination of uniform measures supported on simplices with vertices in S. 



1. Introduction 

The initial motivation for the present paper came from proposed in [16_ efficient 
algorithm recovering an arbitrary convex polytope from axial moments of a poly- 
nomial measure supported on it. This algorithm is based on the formulas for the 
axial moments of polytopes found over 20 years ago independently by M. Brion, 
J. Lawrence, A. Khovanskii-A. Pukhlikov, and A. Barvinok [11|III|IH1[13| , see [HIS] 
for accessible explanation. In [TB] the authors made an essential, although implicit, 
use of a univariate rational generating function for appropriately normalized axial 
moments. Here a multivariate, and explicit, analog of the latter function is devel- 
oped. It turns out it provides a very convenient encoding of non-convex polytopes, 
which is of independent interest. E.g. it leads to a natural definition of vertices 
of such non-convex polytopes, which have similar properties to vertices of convex 
polytopes. It also allows to find the exact solutions of a class of inverse moment 
problems on non-convex polytopes. 

Notation. In what follows we shall always assume that is endowed with a 
fixed coordinate system (xi, a;^), orthonormal with respect to the standard scalar 
product (•, •). Let /i be a finite complex-valued Borel measure in R''. (For standard 
measure-theoretic notions we follow [25\.) Given a multiindex / = (ii, . . . , i^), let 
be the shorthand of the monomial x\'^ . . . x]f and |/| the shorthand for + - ■ ■+id- 
For any multiindex /, define the moment mj (/i) of /i as 

mi{^) := I x^x!^ ...x]^dfi(xi,X2, ■ ■ ■ ,Xd) = / x^(i^(x). (1.1) 

Define the normalized moment generating function -F)j(u) = Ff^(ui, . . . , Ud) of /i 



by 



F^(u):= V Ml^rni iii)u' , where = <... u^^ (1.2) 

ll' ■ • -Id- 

l: = in,...,li)>0 
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Note that F^(u) admits the integral representation 

which is a special case of a Fantappie transformation. For details on the latter, see 
e.g. [U Chapter 3]. A proof of p.3p will be given at the end of Section [5J see also 
Remark [TUl 

Given any complex-valued finite measure fi and any degree 6 homogeneous d- 
variate polynomial p, it is convenient to define the (re)normalized moment generat- 
ing function FP(u) for the measure p/i, where by definition, J^^ fd(pii) = J^j fpdn, 
in such a way that it can be obtained from Ffj_{u) by application of the differential 
operator p(^). Namely, set 

F;(u):^ E M±i^^.Mu^. (1.4) 

l: = ih.....ia)>0 ^' 

Note that FP{u.) ^ Fp^(u) for non-constant p. However, they are also connected, 
by an explicit differential operator as follows. 

Theorem 1. For any complex- valued finite measure p and any homogeneous poly- 
nomial p of degree 6, 

d+S-l / „ 

w= n E-^-a^+^)°^PM(") (1-5) 

£=d \ k ^ 

oF^{u) (1.6) 




p(x)d^(x ) 
(l-(x,u))" 



id+sy.i ,^ '':::zis., . (1.7) 



Here and in what follows o denotes the application of a differential operator to 
a function. The proof of the latter result is basically an exercise in manipulating 
formal power series, and we do not claim its novelty. For the sake of completeness, 
we include a proof in Section [5] 

Results on convex polytopes. A finite set 5 C M"^ is called spanning if it is 
not contained in any (affine) hyperplane in R'^. (Obviously, card(S') > d-l- 1.) As 
usual, by a (compact, convex) poly tope V cM.'^ we mean the convex hull of a finite 
spanning set in W^. The set of vertices of a convex polytope V is the inclusion- 
minimal finite set with convex hull V. A d-simplex in R'* is the convex hull of a 
spanning {d+ l)-tuple of points. By an open polytope (resp. simplex) we mean the 
set of interior points of a compact polytope (resp. simplex). 

Given a convex polytope "P let V = (vi,...,V7v) denote the set of its vertices. 
Assume that V is simple, i.e. each v e V has exactly d incident edges vvg^ , 
vVe_j. Set u'fc(v) := Ve^. — V, for 1 < fc < d. The non-negative real span 
of wi(v),. . . , Wd{'v) is called the tangent cone of V at v. For each define 
detii'vl = I det(wi(v), . . . ,Wd{'v))\ to be the volume of the parallelepiped formed 
by wi{-v),...,Wd{^). 

Given a bounded domain f2 C R'^, we call the measure 

Pn = Xndxidx2 . ..dxd, 

where xn is the characteristic function of fi, the standard measure of fl. 

For a simple convex polytope V we have the following explicit representation of 
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Theorem 2. For an arbitrary simple convex polytope 7^, 

n(,„,(v),u) * • ' 



Remark 1. Instead of the exphcit choice of Wfe(v) for v £ V made above, we can 
take any fixed set of non-zero vectors wi(y)^ . . . , ^^(v), spanning the tangent cone 
of V in P. This does not affect the vahdity of (HH]) and 

Theorem [2] imphes 

Corollary 3. Let A = conv(V) C K*^ be an arbitrary d-simplex. Then 

, , (i!Vol(A) 

n(i-(v,u)) - (i-^^) 

vev 

Remark 2. As we discovered after we proved the above results, statements sim- 
ilar to Corollary [3] in the complex setting can be found in [H Section 3.5] and in 
particular [H Corollary 3.5.6]. 

A variation of (I1.10|) also appears in [5], in the context of designing an efficient 
procedure for integration of polynomials over simplices. 

Notice that an arbitrary convex polytope V admits a triangulation which only 
uses the existing vertices of V, see e.g. [8, Theorem 3.1]. Applying Corollary [3] and 
Theorem [T] to the sum of measures corresponding to such a triangulation we get 
the following. 

Corollary 4. The normalized moment generating function F^{u) of any convex 
polytope V with respect to any homogeneous polynomial density function p of 
degree 5 is a rational function with denominator dividing 

vev 

Example 1. Let A be a triangle in M.'^ with vertices vi = (1, 1), V2 = (2,5) and 
V3 = (3,2). Its normalized moment generating function equals 

7 

Fa{ui,U2) = 7 TT TT r. 

(1 — Ui — U2)[l- — 2ui — 5u2)(l — 3mi — 2U2) 

Its Taylor expansion about the origin up to the terms of degree 7 is given by 

7+42tii+56u2 + 175uj+455ttiU2+329u2+630jii + 2387«i«2-l-3367«i«2 + 1750«2 + 2107M* + 10318ujM2 
+ 21217uiM2 -I- 21546tiiti2 + 8967ti2 + 6762«i + 40082miM2 -I- 106526miM2 + 157976Miti2 + 128772niti2 
+45276u2+21175ui + 145845uiU2+468895«iM2+900123Mitt2-l-10744451«i«2+741993tiiti2+227269ti2, 

which implies that 

7 28 175 455 329 63 2387 

moo = -, "iio = 7, moi = — , m2o = -— , mii = — — , mo2 = — — , mao = — , m2i = — — , 
2 ^ 12 24 12 2 dU 

3591 175 2107 5159 21217 3591 2989 

mi2 = , mo3 = -rriio = , ni3l = , "122 = , mi3 = , "i04 = , 

20 2 30 60 180 20 10 

2863 7609 5642 3066 3025 

"150 = 161, ni4i = , "132 = , '"23 = , mi4 = , mo5 = 1078, meo = , 

15 30 15 5 8 

6945 13397 128589 153493 35333 32467 

msi = ,"142 = ,"133 = ,"124 = ,"115 = , ™06 = ■ 

16 24 160 120 16 8 
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Results on non-convex polytopes. Our second group of results addresses the 
problem of distinguishing different polytopes with the same underlying set of ver- 
tices from information on their moments. The problem of restoring the vertices 
of a polygon or a polytope with a constant mass density from information on its 
moments was addressed earlier in e.g. P^[T5Hr7[|20[|27j . However, the latter do not 
provide the recovery of the vertices in the generality required in the present paper. 
Below we concentrate on the case of constant density and known vertices, and plan 
to return to the general inverse problem for polytopes with unknown polynomial 
density and unknown location of their vertices in the future. 

First we need to define what we mean by a polytope. It turned out that there 
is no general consensus about this notion. Instead there exist several competing 
definitions having their own advantages in different situations. We shall study the 
following class of polytopal objects. 

Definition 1. A subset V CM.'^ coinciding with a finite union of arbitrary convex 
d-dimensional polytopes is called a generalized polytope. 

Definition 2. The number of components of a generalized polytope V is the number 
of connected components of the set V" C V of interior points of V. The closure of 
each connected component of V° is called a component ofV. A generalized polytope 
with one component is called indecomposable. 

Remark 3. We say that a simplicial complex in M"* is pure if all its maximal sini- 
plices have dimension d. Clearly any generalized polytope in can be represented 
as the topological space of an appropriate pure simplicial complex. 

Remark 4. Often one considers a more restricted class of objects, namely polytopes. 
A polytope C M'^ is a generalized polytope homeomorphic to a c?-dimensional 
manifold with boundary. 

We need to introduce the notion of a vertex of a generalized polytope. 

Definition 3. Given a generalized polytope V C M.'^ we call a finite collection of 
open disjoint d-dimensional simplices in R"^ a dissection of V if the closure of their 
union coincides with V. 

A wealth of material on dissections of polytopes can be found in j21] , see also 

Definition 4. Given a generalized polytope 7-* C K'' we call a point v a vertex of 
V, if V is a vertex of (the closure of) some open simplex in every dissection of V. 

Definition 5. Given a point p € V of a. generalized polytope V we denote by the 
tangent cone Tp{V) of 7-" at p the set obtained as follows. For a sufficiently small 
e > set 'Pp(e) = "P n Bp{e) where Bp{e) is the e-ball centered at p. Define Tp{P) 
as the set obtained by taking a ray though p and every point of Vp{e). In other 
words, TpiV) is the cone with the apex at p and the base Bp{e). (Obviously, Tp{V) 
is independent of e for a sufficiently small e > 0.) 

Lemma 5. A point v is a vertex of V if and only if T^{V) does not admit a 
decomposition in the disjoint union of convex polygonal subcones, such that each 
subcone in the decomposition has a translation-invariant direction. In particular, if 
the tangent cone to P at v has a connected component with no translation-invariant 
direction then v is a vertex. 

We denote by conv(S') the convex hull of an arbitrary set S C W^. The above 
lemma implies that any vertex of conv(P) is a vertex of V. 

The following result extends Corollary [3] to the case of generalized polytopes. 
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Proposition 6. For any generalized polytope V with the set of vertices V('P), the 
denominator of its normahzed moment generating function Fy,{\i) with respect to 
a homogeneous polynomial density function p of degree 5 divides 

<I>p(u):= n (l-(v,u)r'. 
vevCP) 

Remark 5. There exist generalized polytopes which do not admit dissections with 
only existing vertices. The simplest example of this kind is the Schonhardt polyhe- 
dron, see Figure [T] and [5^. Absence of a dissection T which uses only its 6 vertices 
can be established by observing that none of the edges AC, A'B, and B'C can 
appear in a simplex of T, yet any simplex on these 6 vertices must contain one of 
them. Therefore, Proposition IH] is not an immediate consequence of Corollary [31 




Figure 1. Schonhardt polyhedron obtained from an octahedron 
(on the left) by removing tetrahedra [ABB'C], [AA'B'C], and 
[A'BCC"]. 

Remark 6. For "generic" generalized polytopes V the denominator ri(u) of i^-p(u) 
equals $75 (u), but for certain special polytopes the denominator fl(u) may be its 
proper divisor, as can be seen from the following example. Let A — {0, ai, 02, 03} C 
be a spanning set, and w e M^. Let 'P± := conv(w ± A) and V := V+UV-. 
Then 1 — {u,v) does not appear in il{u), as 

E (u,a,:)(u,aj) + (1 - (u,?;))2 

Fp(u) = F-p^ (u) + Fr_ (u) - A'i^^^^^i— — 

l<i<3 

where iiT ^ is a real constant. 

Now we introduce several finite-dimensional linear spaces related to a given finite 
spanning set 5* C K'^. Let V{S) be the set of all generalized polytopes V whose sets 
V{'P) of vertices are contained in S. For V G 7^(5') we denote by /z-p its standard 
measure. (Obviously, fi-p is supported on "P C conv(S').) 

Denote by 9Jt(S') the linear space of all signed measures, i.e. the linear span of all 
standard measures /Lt-p for V S ViS). Let 9Jt^(S') C 931(5') be its subspace spanned 
by Ma, for A G 'P(S') a d-dimensional simplex. (The space 971^(5) has earlier 
appeared in [T] , [2] , [3] in a somewhat different context.) We shall refer to elements 
of 9)t(S') as to polytopal measures with the vertex set S. The next conjecture is 
central to our study. 



6 



N. GRAVIN, D.V. PASECHNIK, B. SHAPIRO, AND M. SHAPIRO 



Conjecture 7. For an arbitrary spanning set S and any generalized polytope V 
with set of vertices contained in S, its standard measure /ip belongs to m^{S). In 
other words, m{S) = m^{S). 

By Remark [5l the above conjecture is non-trivial. While we do not have a 
proof of Conjecture [7] in its full generality, we have succeeded in proving it for a 
rather large class of spanning sets. Roughly speaking, the latter should be close to 
"generic". Specifically, given a finite spanning set 5 C K.'^, we say that S is weakly 
non-degenerate if any (d + 2)-tuple of points from S is spanning. If S satisfies the 
stronger condition that each {d + l)-subset of S is spanning then we call the latter 
S strongly non-degenerate. 

Theorem 8. Conjecture [7] holds for any weakly non-degenerate finite set S. 

Remark 7. Theorem [8] would imply Conjecture [7] if one could prove that the 
standard measure of an arbitrary generalized polytope V can be obtained as the 
limit of the standard measures of a 1-parameter family of generalized polytopes Vit) 
with V{0) — V such that for t ^ the vertices of V{t) are weakly non-degenerate, 
and each vertex of V{t) tending to a vertex of 'P as t 0. Unfortunately we are 
unable to prove the existence of such deformations in general. 

The key idea in the proof of Theorem [S] is to study the corresponding spaces of 
Fantappie transformations of signed measures in 9Jt(>5'). In particular, we are able 
to compute the corresponding dimensiontQ. In more detail, let ^{S) (resp. "^^{S)) 
be the linear space of Fantappie transformations of signed measures in 371(5') (resp. 
dR^{S)). In other words, 'S{S) (resp. ^^{S)) is the space of normalized moment 
generating functions of signed measures in 971(5) (resp. 9Jl'^(5')). 

Since each compactly supported measure is uniquely determined by its complete 
set of moments, the map 

i^^ : 971(5) ^ y(5), (1.11) 

induced by the Fantappie transformation is a linear isomorphism, cf. |4i Section 
3.5]. 

Finally, given a spanning set 5 = {vi,...,vjv} C M'', denote by 9^at(5) the 
linear space of all rational functions with the denominator $s(u) as in (|1.12[) . 

N 

$5(u)=n(l-(v.,u)), (1.12) 

1=1 

and with the numerator an arbitrary real (inhomogeneous) polynomial of degree 
at most N — d — 1. Here the numerator and the denominator might have common 
factors. 

Proposition 9. ^^{S) coincides with 9^at(5) if and only if S is strongly non- 
degenerate. 

Corollary 10. If 5 is strongly non- degenerate then 97l'^(5) = 971(5). 

Corollary [10] implies that for strongly non-degenerate 5 the dimension of all 
these linear spaces equals (^J^)- Note that Corollary [Till settles Theorem [5] for the 
strongly non-degenerate 5. 

Our final goal is to explicitly solve the following inverse moment problem. 



Note that presently we are not aware of a formula or a recipe for calculating the dimension 
of 2K^(S') without the assumption of the theorem. (The manuscript [2] contains an algorithm 
constructing a basis of this space.) 
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Problem 1. Given a strongly non-degenerate spanning set S* C M"*, jS*] = A^, find 
the unique polytopal measure in 971(5') with a given set of all moments up to order 
N -d-l. 

We start with the following simple observation. 

Lemma 11. Given an arbitrary spanning set S C R'^, \S\ = N, and an arbitrary 
polynomial r(u) of degree at most N — d—1 there exists a unique rational function 
i?(u) = P(u)/$5(u) with Taylor polynomial of degree N — d — 1 at the origin 
equal to r(u). Namely, P{u) = [T(u)$5(u)]^_^_-^, where [-Ijv-d-i stands for the 
truncated polynomial with all monomials up to degree N — d — 1. 

For S = {vi, . . . , vjv} C strongly non-degenerate, we give an explicit inver- 
sion formula determining the densities of an unknown polytopal measure having a 
given set of moments up to order iV — d — 1 on each simplex in a natural basis of 
9Jt'^ [S). In view of Lemma lTTI we can assume that we are already given an arbitrary 
rational function R{u) = P(u)/$5(u), where degP(u) < N — d — 1, and we want 
to determine the densities of the required signed measure from 971(5) in terms of 
numerator -P(u). 

From now on we shall choose the basis of 97t'^(5) consisting of the standard 
measures of all simplices containing the last vertex v^r, see Lemma Il5l below. Let 
•2 = {^1, h, ^jv-i} be the {N — l)-tuple of linear forms corresponding to vertices 
Vi, V2, . . . , vjv-i, where li{u.) = 1 — (v,i,u). Consider the linear span Vo of all 
possible products of the form Ij-^ ■ Ij^ ■ ... ■ 1 < ji < j2 < ■•• < Jjv-d-i- 

There are (^^^) such products, and each of them is a polynomial of degree at most 
N — d — 1. On the other hand, the dimension of the space Pol{N — d — l,d) of 
all (inhomogeneous) polynomials of degree at most N — d — 1 in d variables equals 
(^-^), as well. 

Define the square matrix Mats of size (^J^) with entries being coefficients of 
the above products of linear forms with respect to the standard monomial basis in 
Pol{N — d — 1, d). We assume that Mats acts on the space Vs, of column vectors. 

Theorem 12. For an arbitrary strongly non-degenerate spanning set S C M'', 
\S\ — N, the matrix Mats is invertible. Moreover, for a rational function i?(u) = 
P(u)/<I'5(u), where P(u) is an arbitrary polynomial of degree N — d — 1, there 
exists a unique measure fiij G 971(5) with Fantappie transform i?(u). Namely, 

fiR = Matg^{P{u)). (1.13) 

Remark 8. A detailed explanation of the meaning of (|1.13p can be found in the 
proof of Theorem II 2[ see also Example [5] below. An explicit formula for the matrix 
Alatg^ is given in Lemma [TSl 

Recall that a spanning set 5 is weakly non-degenerate if any {d + 2)-tuple of its 
points is spanning. With minor changes, the above solution of the inverse moment 
problem can be adapted to this more general case. In order not to overload the 
introduction we refer the readers interested in this situation to Section |31 The 
case of an arbitrary spanning set 5, however, remains unsolved and offers several 
interesting challenges in matroid theory. We hope to return to it in the future. 

It will be convenient to work with scaled volumes of simplices, which we call 
weights. 

Definition 6. Given a signed measure /i in R'' and a d-dimensional simplex A C M'' 
we define the weight wa of A by the formula: 

WA=dl / dfi. (1.14) 
J A 
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In other words, the density c?a of the measure in question which should be placed 
at A equals 

^ " d!Vol(A) • 

We finish the introduction by explicitly solving the above inverse problem for a 
concrete 5-tuple of points in M^. 

Example 2. Set S — {vi, V2, V3, V4, V5} where vi = (1,0), V2 = (2,1),V3 = 
(1,2),V4 = (0,1), V5 — (0,0). The corresponding set £ = {^i, ^2, ^3, '4} of linear 
forms is given hy h — 1 ~ ui, I2 — 1 — 2ui — M2, I3 — 1 — ui — 2u2, — 1 — U2- 
Additionally, = 1. We are considering the basis of SR^{S) consisting of (the 
standard measures of) 6 triangles containing V5. Therefore we need 6 quadratic 
forms obtained as pairwise products lilj, 1 < « < j < 4. We get 



'hh = 


1 


— 3ui 


- M2 H 


h 2u\ + U1U2 




1 


- 2ui 


- 2m2 


+ u\ + 2uiU2 


hh — 


1 


-Ul~ 


' U2 + 


UIU2 


hh^ 


1 


— 3ui 


- 3U2 


+ 2u\ + 5uiU2 + 2u2 


hh = 


1 


- 2ui 


-2U2 


+ 2uiM2 + U2 


^3^4 = 


1 


- Ul - 


- 3u2 H 


h U1U2 + 2u\. 



Notice that hh corresponds to triangle A345, hh to A245, hh to A245, hh to 
A234, hh to A145, hh to A135, and hh to A125. Ordering monomials spanning 
the space Pol{2, 2) as (1, ui, U2, uf, U1U2, u\) we get the 6 x 6-matrix Mats and its 
inverse AlatZ^ as follows 









Mat 


s = 












AMat 


-1 

s — 








hh 


hh 


hh 


hh 


hh 


hh 




1 


Ul 


U2 


ul 


U1U2 


ul 


1 


1 1 


1 


1 


1 


1 


1\ 


hh 




-1 


1 


1 


1 




Ul 


-3 


-2 


-1 


-3 


-2 


-1 


hh 







-4 








-4 


U2 


-1 


-2 


-1 


-3 


-2 


-3 


hh 


9 


3 


3 


1 


1 


1 


ul 


2 


1 





2 








hh 


1 


1 


1 


1 


1 


1 


UIU2 


1 


2 


1 


5 


2 


1 


hh 


-4 


-4 





-4 








ul 


I 








2 


1 


2) 


hh 


I 1 


1 


-1 


1 


-1 


iJ 



(For T^nical reasons we give iAIatg^ above.) Thus, given an arbitrary rational 
function i?(ui, M2) — P(ui, U2)/$s(wi, U2) where P(ui, U2) = aoo + 01^0^1 + ao, 1^2 + 
a2.o'^i + aiiUiU2 + ao2ul is a polynomial of degree at most 2 and ^s{ui,U2) = 
hhhhh, we get 



W345 = 


i:(aoo - 


- flio + floi + a2o 


- On 


+ 002) 


W245 = 










W235 = 


i(9aoo 


+ 3aoi + 3aio + 


020 + 


On + O02) 


Wl45 = 


j(aoo -f 


- aoi + aio + 020 


+ On 


+ 002) 


W135 = 


—aoo — 


OlO ^ 020 






.W125 = 


j{aoo -f 


- aio — aoi + 020 


- On 


+ 002), 



where Wijk is the weight of the signed measure to be placed on Aijk , see (|1.14p . 

To illustrate all steps of solution of our inverse moment problem assume that 
we are looking for a polygonal measure with the vertex set S and (ad hoc chosen) 
moments moo = 1,^10 = 2,moi = 3, 77120 = 4, mn — 5,mo2 = 6. Then its 
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normalized moment generating function -F)i(u) satisfies tlie relation 



F^^in) = 1 



2! 



3! 



3! 



4! 



4! 



4! 



P{ui,U2) 



0!0! 1!0! 0!f! 2!0! ' f!l! 0!2! ^ hhhUk 

where P{ui, U2) is a (non-homogeneous) polynomial of at most second degree. Thus, 
truncating the product of the left-hand side and hhhlih up to the second degree, 
we obtain 

P{ui,U2) = 2 + 4ui + 10u2 + lOul + 24uiU2 + IOU2, 
i.e. floo = 2,aio = 4, oqi = 10, 020 = 10, an = 24, ao2 = 10. Thus ^345 = 
1,^245 = — 22, W235 = 26,i(Ji45 — 15,u'i35 = —16,^125 = "2. The areas of the cor- 
responding triangles are equal to: Area(A345) — i; Area(A245) = 1; ylrea(A235) = 
|;Area(Ai45) = i; ylrea(Ai35) — l;^rea(Ai25) = ^. This implies that the den- 
sities of the measure of the corresponding triangles are equal to c?345 — 1 , ^245 = 
— 11,^235 — ^,di45 — 15,(ii35 — — 8, di25 = —2. To obtain the final densities in 
the convex hull conv(S') of S one has to decompose conv(iS') into domains obtained 
by removing from conv(5') the set of all hyperplanes spanned by vertices in S. For 
each such domain we should add up the densities of all basic simplices containing 
this domain. The resulting measure is shown in Fig. [51 



V4-(0,l) 




V2 = (2,l) 



V5 = (0,0) 



vi = (l,0) 



Figure 2. Final measure in Example 



Remark 9. Domains into which the convex hull conv(5') is cut by the hyperplanes 
spanned by S were introduced in [T] where they were called chambers. The incidence 
matrix of the simplices spanned by S and those chambers was studied in some 
detail in [2], [3]. This matrix allows to formalize the last step of construction of 
the above polygonal measure, where information on the densities of the simplices 
is transformed into information on the densities of the chambers. But, in general, 
already the number of chambers is a complicated invariant of the set 5*. It seems 
that the general problem of constructing the set of chambers and the corresponding 
incidence matrix in terms of a given S is quite non-trivial. 
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2. Proving results on convex polytopes 

Following Brion-Lawrence-Khovanskii-Pukhlikov-Barvinok, see [HIHlIinillH] i we 

define for each vector z S K'^ the j-th axial moment /ij(z) of a simple convex 
polytope V with respect to z as 

lij{z) = / (x, z)^ (ix. 
Jv 

We will use the following important statement, cf. e.g. [5^, Theorem 10.5]. 
Theorem 13. The moment /ij(z) satisfies 

= mTTE^'^)'^'^v(z), (2.1) 

^■^ ^' vGV 

where D^lz) I dctK^\ — ^ arbitrary vector for which the products 

^ ' nj=i(«'j(v),z) ' ^ 

n^=i('"^i('^)i^)' V 6 V, do not vanish. Moreover, the following identities hold: 

^(v,z)^Dv(z) = 0, forO<j<(i-l. (2.2) 
vev 

Proof of Theorem [H To prove p.8p , consider the generating function 

ij + dy. 



where m £ R. Formula (|2.ip implies that (u) is rational. Indeed, 

j=o vev n/c=iK(v),z) 
-( D'Y V (v,u)'^|detj^.| (-1)-^ 

v'bntiK(v),z)^eo ' v^nLiK(v),u) l-(v,u)' 

where u = uz. On the other hand, using the multinomial coefficients ('^') = ji/^'j^] 
of multiindices J = (ji, . . . , jd) l~ | J|, one gets 
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where mj['P) = m^di-p). Therefore, 



F-p{uz) := F-p{uzi,...,uzd) ■■=^ ^ 



E 

3=0 



E 

•/: = (jl,... Jd)l-j 



Jl! • ■ • Jd! 



and follows. 

In view of relations (|2.2p the right-hand side of (|1.8p can be rewritten as (|1.9p . 
Indeed, writing (1 — (v, u))~^ = Ejlo('^'^)"' ^'^'^ expanding (|1.8p with respect to 
jth powers of (v, u), we see that (|2.2p implies that for j < d the sum of all terms 
(v, u)-' vanishes. □ 



Proof of Corollaryl^ Let V = (vq, vi, . . . , v^). Then for each j i we have 
Wj(vi) = Vj — Vi. Hence | detifvj does not depend upon i and equals (i!Vol(A). 
The right-hand side of (jl.Sp becomes 



^lfd\Yo\{A)y —J^^ i , 



= (-l)'^d!Vol(A)^- 



(v,,u)" 



^ nj5:^»(vj - v,,u) l-(v,,u) 
where = (v^, u). Computing the common denominator of the latter, we obtain 



i^A(u) 



d 


n (a -con (1-0) 




(-l)''rf!Vol(A) ^=o 




nto(i-c.) 


Us>tiCs-Ct) 





It is convenient to introduce one more linear form Cd+i 1, so that the last 
expression reads as 



Fa{u) = (-l)''d!Vol(A)- 



d 

E 


n (a - Ci) 






d+l>k>l>0, 











n (c. 

d+l>s>t>0 



CO 



(2.3) 
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To complete the proof we notice that 



= det 



/I 


1 . 


. 1 


i\ 


Co 


Ci ■ 


• Cd 


1 


Co^ 


a ■ 


■ a 


1 


c^' 


Ct ■ 


■ ci 


1 


VCo' 


Cf ■ 


■ c.^ 


1/ 



(-l)'+'Ecf(-ir-det 



i=0 





/I 1 ... 








Co Ci ■■■ 






1 • det 


Co^ Ci^ ... 


a 


+ 




Uo' cf ... 


a) 




/I .. 


1 1 






Co . 


• Ci-1 Ci+1 




Cd+l 


Co^ . 

1.^ 


. 'ii-l 

^.1 L..J 




^.1 



, /-d /-d /-d /-d 



= n (Cfe-o)+(-i)'+'Ec'(-i)' 



d>fc>/>0 



d+l>k>l>0. 



Indeed, the first matrix has two identical rows and thus vanishing determinant, 
which we expand with respect to the last row. The last equality is the standard 
formula for the Vandermonde determinant. Thus we have 



i=0 



n (Cfc-o: 



d+l>fc>/>0, 



d>k>l>0 



Now we plug this formula into (|2.3p and get 

1 



Fa(u) = d!Vol(A)- 



n(i-o.) 

t=o 



d!VoI(A) 

n(i-(v,u))- 

vev 



□ 



Lemma 14. Let u — {ui, . . . ,Ud) and x = (xi, . . . ,Xd) be formal variables, and 
eeR . Then 

(E + ^) ° (1 - = - (2.4) 

Proo/. Note that Ufegf;: o (1 - (x, u))-^ = XkUk£{l - (x, u))-^-i. Thus 

(E + ° (1 - (x, u))-^ = ^(1 - (x, u))-^-i(x, u) + £{1 - (x, u))-^ 



£(l-(x,u))- 



□ 



Proof of (|1.3p . For a d-variate polynomial g{z) we denote by g (u^) the differen- 
tial operator g (^ui-^, ■ ■ ■ , ""dgf^^ . We use the identity 



9 I ) ° E "I'^^u^ = E 



(2.5) 



which holds for any formal d-variate power series a/x^u^ and any c?-variate 
polynomial g{z). (It can be easily verified for monomial g(z) and then extended by 
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linearity.) Setting /i(z) := HLi (ELi Zk + ij , notice that h{I) = (|/| + 1)(|/ 
2) • ■ • (|/| + d). Now using (12. 5p together with the obvious identity: 

(l-(x,u)r=^(l^l)xV, 

one obtains 



/>0 



/l u 



(i/i(x) 
9uy 1 — (x, u) 



dl dfi{x.) 



(l-(x,u))'^+i' 



where in the final derivation we repeatedly made use of (12. 4p , for 1 < £ < d. □ 

Remark 10. Another point of view on (II. 3p is that it is the result of the application 
of the differential operator g (u^) to the integral transformation J^^ il^i^u) 
measure /i (also known as the Fantappie transform of /i); see e.g. |19j . 

In |22J a similar idea was applied to the harmonic polygonal measures in the 
plane. 



Proof of Theorem\^ Assume first that p(a;i, = x 



K 



■ X, IS a mono- 



mial and consider 



m + d)\ 



E 



m/(^)u^. 



One gets 



i\K\ 



( 



\ 



K 



I={ii,...Aa)>a 



{\ I\ + \K\+d)\ ^{i,+k,)\ 

= E 

/=(ii,...,id)>0 



J 



ail*.! 



Observe that the normalizing coefficients of mj{x^ /j.) in the latter expression de- 
pend only on / and \K\ but not on particular entries of K. Therefore for an 
arbitrary homogeneous p of degree S one gets by additivity 



p(|i)o^.(u) 



E 

/=(n,...,id)>o 



i\I\+d + S)l 



d 

n 



This shows (|1.6p . Repeated application of (|2.4p . for d + 1 < £ < d + 6, to the 
integral representation (|1.3p . respectively, to the representation p.2p . of i^pp(u) 
implies (|1.7p . respectively, p.5p . □ 
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3. Inverse moment problem for strongly non-degenerate S 

Proof of Lemma [31 We prove first that the tangent cone at any non- vertex aUows 
a decomposition into convex polytopal cones each having a translation-invariant 
direction. 

Let V be a point in V which is not a vertex. Then there is a dissection T of 

V such that v is not a vertex of any simplex of T- Let U be the set of simplices 
Su of T with closures containing v. Take the dissection of the tangent cone Ty{V) 
into the tangent cones to simplices from [/, T^{V) — Uu£uTv{Su)- Clearly, every 
subcone r„(5„) contains a translation-invariant direction (any direction parallel to 
the minimal face containing v). 

Vice versa, to prove the converse implication, let us take a dissection of the 
tangent cone T^{'P) into a disjoint union of convex polytopal cones Qi, . . . , Qk- By 
definition of the tangent cone and since V can be represented as a finite union of 
simplices we obtain that any sufficiently small neighborhood of v in the tangent 
cone T^{V) is a neighborhood of v in the entire V. Consider the parallelepiped 
Boxg centered at v that is e ball centered at v, in the ii-norm. Note that each 
convex polytopal set Qi n Boxg can be decomposed into a union of simplices that 
do not contain v as a vertex. 

Further notice that the set V \ Box^ can be represented as a finite disjoint union 
of simplices, since Box^ is the intersection of a finite number of half-spaces and V is 
a disjoint union of simplices. Clearly, every simplex in this union should not have 

V as a vertex. Now combining the dissections of each Qi fl Box^ and V \ Box^ we 
obtain the required dissection of T-". □ 

Proof of Proposition We begin by considering the case p = 1 . Let 7" be a dissec- 
tion of V with vertices V(7~). Corollary [3] implies that ^^(u) has a denominator 
dividing 5r(u) = nvev(r)(l - (v,u)). Take Vi G V(r) \ V{V). Then there ex- 
ists another dissection T' such that vi ^ V{T'). Expressing F-p{u) as ratios of 
polynomials, we have 

Here g-r' is not divisible by 1 — (vi, u), by the choice of 7"'. Thus fj- is divisible by 
1 — (vi,u), and can be canceled out in the expression for F-p{\i). 

The case of arbitrary homogeneous p follows immediately by applying Theorem[T] 
to the already covered case p = 1. □ 

Proof of Proposition First we show that for an arbitrary finite spanning set 
S" C M'' the space 9Jl^ (5*) has a basis of d-dimensional simplices containing a fixed 
vertex v € S. In particular, the set of all d-dimensional simplices containing v 
spans 9)1'^ {S) but is not necessarily a basis. Consequently, their Fantappie trans- 
formations spans 'S'^{S). The following result is formulated as Theorem 4.2 of [5] 
and in a different form in [Ij. (We omit the proof of this statement here.) 

Given two points p and q and a set M in W^, we say that q is visible from p with 
respect to M if the straight interval pq is disjoint from M. 

Lemma 15. Given a d-dimensional simplex a C M.'^ denote by V{a) the set of 
vertices of a. Let cr" be the interior of a. Let p be any point in K'^ and let 
(resp. Q~) be the set of all (cZ— l)-dimensional faces of a which are visible (resp. not 
visible) from p with respect to (T°. Then the standard measures of all d-dimensional 
simplices with vertices in V(cr) U {p} satisfy 
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where /io-.^p is the standard measure of the d-dimensional simplex spanned by the 
vertices of cji and the point p. 

Remark 11. If degenerate simplex, i.e. p lies in the hyperplane spanned 

by CTi, we simply exclude the corresponding term ^ from the above formula. 

To prove Proposition [51 we need to show that '^^{S) coincides with 5Hat(iS') if 
and only if S is strongly non-degenerate. Indeed, 5'^(S') C 9^01(5*) for an arbitrary 
spanning S, by Proposition [SI The Fantappie transform : 971'^ (S*) S^iS) is 
a linear isomorphism which implies that dim 9}l^(S') = dim5^^(S'). By Lemma [T51 
the space 371'^ (5*) is spanned by the standard measures /ia of the set Bi of all d- 
dimensional simplices containing the fixed vertex v^. Let us fix the vertex v„ and 
consider the set For S strongly non-degenerate the cardinality of i3„ equals 

Now we show that d'^iS) — 5Hat(5), where 9iat(5) has the dimension (^^ as 
it is isomorphic to the space Pol{N — d—l,d) of all d-variate polynomials of degree 
at most N — d — 1 . This would immediately imply that the standard measures of 
simplices in y8„ are linearly independent. 

Lemma 16. If 5* is strongly non-degenerate, then 'S^{S) — D\ai{S). 

Proof. WerecaUthat5-'^(S') comprises all linear combinations of the rational func- 
tions 

rf!Vol(v,;i,...,v»^,v») 

(1-K,u)) (l-(v,,,u))-(l-(v„,u))- 

For each term 1 — (vi,u) we consider a (homogeneous) linear form li{uo,u) = 
uq — (vi, u) in d-l- 1 variables uq, . . . , Ud- Set n = N — 1 where n > d+1. For the n- 
tuple £ = {^1, 12, In} of linear (d-f l)-variate forms, let be the linear span of all 
possible products of the form h-^li^ ■ ■ ■ ^i„_d, 1 < ii < 12 < •■• < in-d < n. Observe 
that Vc is the space of all numerators that one can obtain in ^'^{S). We need to 
show that Vfi contains HPoly{n~d, d+1), the space of all d+l-variate homogeneous 
polynomials of degree n — d. Recall that any d+ 1-tuple of linear forms , . . . , h^^^ 
is linearly independent due to the strong degeneracy assumption. Thus we can 
express each single variable uq, . . . , as a linear combination of these forms. Since 
Vc contains all products kji^ ■ ■ where j G {1, ...,n} \ {ii, . . . ,i„_d_i}, 

we conclude that contains all homogeneous polynomials of the form 

hji2 ■ ■ ■h„_a-i'^k, for < fc < d. 

From that we deduce that V£ contains all homogeneous polynomials of the form 
hj^k^ ■ ■ ■ h„_a-2'^kUj, where j, k € [n] \ {ii, . . . , in-d-i}- Continuing along the same 
lines we derive by induction that Vc contains HPoly{n — d,d + 1). □ 

For an arbitrary spanning S the cardinality of i3„ is at most (^J^) = dimlHat(5'). 
Furthermore, if S is not strongly non-degenerate the cardinality of Bj^ is strictly 
smaller than ( J^), as there will be linear dependencies among the standard mea- 
sures on the simplices in B^. Therefore dim3^^(S') < dim9lat(S'). □ 

We define the square matrix Mats, of size (^) with entries being coefficients of 
the above products of linear forms w.r.t. the standard monomial basis in HPol(n — 
d,d+l). 

Lemma 17. The determinant of Mats, is proportional to the product of the de- 
terminants of all (d + l)-tuples {li-^jli^, . . . ii < 12 < ... < id+i- (By the 
determinant oi a {d + l)-tuple of vectors in M.'^^^ with a fixed basis we mean the 
determinant of the matrix formed by the coordinates of these vectors in a chosen 
basis.) 
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Proof. Indeed, det{Mat£,) is a form of degree (d + 1)(^) in the coefficients of the 
finear forms ^i,...,Z„. Thus the product Yii^ ^ det(Zij , i • ■ • j ^id+i) the 
same degree as det (Maio). Therefore it suffices to show that det (Mat o) vanishes 
as soon as some of det(/ij , ^^37 ■ • ■ i 'id+i) vanishes. (Observe that all polynomials 
det(Mai£) are coprime.) Without loss of generality assume that is a linear 
combination of I2, ■ ■ ■ ,ld+i- But then the column of Mat^, corresponding to the 
(n— o?)-tuple (1, d+2, d+'i, . . . , n) will be a linear combination of those corresponding 
to (2,d + 2,d + 3,...,n), ...(d+l,d + 2,d + 3,...,n). □ 



Proof of Corollary] 101 As we mentioned above, dJl{S) is isomorphic to ^{S) and, 
analogously, DJl^{S) is isomorphic to ^^{S). Thus, if we prove the equality ^{S) = 
d^{S) then we get m{S) = m^{S). By Lemma[ini the space ^''^(5) coincides with 
the linear space of all rational functions with the numerator an arbitrary polynomial 
of degree at most N — d — 1 and the denominator <i>5(u) equal to the product of 
all linear forms dual to all vertices in 5*. By Proposition [S] an arbitrary function 
in ^{S) is a rational function with denominator of desired form and numerator of 
degree at most N — d—l, for obvious reasons — take an arbitrary dissection and sum 
over its simplices. Since all such functions are already in d^{S) we are done. □ 



Proof of TheoremM'A Given a strongly non-degenerate set S ~ {vi, . . . , v^r-i, vjv} 
and the Fantappie transform i?(u) ~ P(u)/$5(u), where <i>5(u) — JI^Li^j^u), 
we want to solve the inverse moment problem. (It is easy to obtain P(u) from 
information on the moments of order at most N — d — 1 using Lemma [TT]) 

To solve the latter inverse problem using Corollary [3l we need to find an appro- 
priate set of weights w = ii < ii < . . . < where is the weight 
(recall Definition I1.14p of the d-dimensional simplex conv(vjj^ , v,i2 , . . . , v^^, v^v) so 
that 



E 



^ ^ii(")^»2(u)---'-M(u)^«(u) ;i(u)?2(u).../„(u)' 

Clearing the denominators, we get the equation 

E W^,^...,,Jj,{u)lj^{u) . . .lj^_^_^{u)lN{u) =P(u), 
il<i2<---<id 

where {ji, . . . j'jv-d-i} = {1, 2, . . . , A'^ — 1} \ {ii, . . . , i^}. The latter equation is 
obviously equivalent to the system of linear equations 

Mats ■ w = (p/j, . .■,pit), where P(u) = ^pju' , 

I 

and w = is the vector consisting of the weights of all simplices containing 

Vat. □ 



Theorem[T^ solves the inverse moment problem for strongly non-degenerate span- 
ning set S. We can make this solution more explicit by giving a closed formula for 
the inverse matrix Alatg^. To do this we introduce an extra variable wq G M and 
identify the space Pol{N — d—l,d) with the space HPol{N — d—l,d^\) of homo- 
geneous forms of degree N ~d~l iiid+l variables {uq,ui, . . . ,Ud)- We homogenize 
each linear form li{u) in £ as /i(u,Mo) — uo — (vi,u). (The matrix Mats remains 
unchanged.) 
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We also need the following {d + I) x {N — 1) matrix L 

h h ■ ■ ■ In-1 

I \ 1 ... 1 



L 



-(vi,ei) -(v2,ei) 



Ud \-(vi,erf) -(v2,ed) 



-(vAr_i,ei) 



'(vAr_i,ed) / 



(3.1) 



associated with Matg. For all possible subsets = {ii, . . . , i^}, i\ < < ■ ■ ■ < 
of d distinct columns of L consider the linear in u function Li[^] (u) given by: 



(u) = det 



Wo 



1 

(v,i,ei) 



Ud -(ViijCd) 



-(v»<i,ei} 



(3.2) 



Denote by Lij^jj j the coefficient of Uj in the linear form Li[(i](u). For each 1 < j < 
N - \ define L(j,i[d]) as 



L(j,i[. 



dot 



1 



(3.3) 



Note that if j ^ then L(j, i[(i]) ^ 0, as by the assumption of strong non- 
degeneracy of S the corresponding d + 1 linear forms are linearly independent. On 
the other hand, if j G then we have L(j, i[(i]) = 0. 

The matrix Matg has the following explicit description. 

Lemma 18. For each {N ~ d — l)-tuple of forms {/j^, Ij^, . . . , Ij^-d-A set i[d] = 
{zi, . . . ,id} ^ [N - 1] \ {ji, . . . jN-d-i}, where ii < 12 • • • < id- Then, 

nn ni 

Un Wl ... 



JV-d-1 
k = l 



\ 



(3.4) 



Proof of Lemma [751 In order to show that Matg ^ defined by p. 41) is indeed the 
inverse of Mats we need to verify that Matg^ ■ Mats is the identity operator on 

Let e' be the standard basis vector of corresponding to the product of linear 
forms L> . . . L' . Then Mats ■ e' is the vector consisting of the monomial coefE- 
cients of the homogeneous form Ij'^ . . . Iji^ ^ ^ in the variables uq, . . . ,Ud. Let e-^ be 
the row vector of Alatg^ corresponding to the product Ij-^ . . . ljf^_j_i- We note that 

N-d-l 

in e-^ • {Mats ■ e') one can factor out the common denominator J| L(jfe, i[d]) of 



all fractions in e^; the remaining factor is of the form 



k=l 



Mat5[u^e']nL;;;i],, ='j;(Li[<i](u))...Z,;,_,_,(Li[,](u)). (3.5) 



I=(nQ,...,nd) 



3=0 
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Note that = L(j^,i[d]) for 1 < fc < TV - d - 1, i.e. this holds for 

all terms in the product on the right-hand side of p.Sp . Hence, if e 7^ e' then 
-{M ats ■ e') = 0, as among j'l, . . . , i/v-d-i o"^^ '^^^ ^^^'^ f ^ i[c?] with L(j', \[d]) = 
0. On the other hand, if e and e' coincide, then the right-hand side of (I3.5P is equal 

N-d-l 

to W L(jfc,i[(i]). Dividing by the common denominator of the fractions in e, we 

fe=i 

obtain ■ {Mats • e) = 1. □ 



4. Inverse moment problem for weakly non-degenerate S 

Given an arbitrary spanning set S — {vi, V2, . . . , vat}, consider the linear space 
6(5') C Pol{N ~d~l,d) spanned by all products Ij^ Ij^ ■ ■ ■ 'jjv-d-i ' Ji < J2 < • ■ • < 
jN-d-i- The next statement explains why we can extend our solution of the inverse 
moment problem from the case of strongly non-degenerate S to the case of weakly 
non-degenerate S. 

Lemma 19. 0(5*) ~ Pol{N — d — l,d) if and only if S is weakly non-degenerate, 
i.e. each (d + 2)-tuple of points of S is spanning. 

Proof. We have N (non- homogeneous) linear forms Ii...,In in variables u = 
(ui, . . . ,Uii) and the linear space spanned by all possible products of {N — 
d — l)-tuples of distinct forms. We need to investigate whether coincides with 
Pol{N — d — l,d). Homogenizing, we consider the same question for the linear ho- 
mogeneous forms and the homogeneous polynomials of degree N — d—1 in variables 
(uo,ui,...,Ud). 

First assume that there are d+2 linear forms li, . . . , ld+2 which are not spanning. 
Then one can find a non-zero vector zq £ E'*+^, such that li{z) = ■ ■ ■ — ld+2{z) = 
0. Note that each product of — 0? — 1 different forms chosen from li, ... ,1m 
contains at least one form among {Zi, . . . , ^^+2}- Therefore any linear combination 
of products of — 0? — 1 forms vanishes at Zq. Thus Vs, cannot coincide with 
HPol{N -d-l,d+l). 

Conversely, assume that every (d -I- 2)-tuple of distinct forms among Zi, . . . , /jv 
is spanning. First, we notice that HPol{N — d — 1, d + 1) can be spanned by the 
all possible products of iV — d — 1 linear forms (not necessarily pairwise distinct). 
Indeed, since first d+2 forms span the dual space of M''"'"^, we can express each 
variable Xi cLS Ql, linear combination of these forms. Therefore every monomial of 
degree N — d—1 can be expressed as a linear combination of products of iV — d — 1 
forms. 

Now we show that each product of — d — 1, not necessarily distinct, forms 
can be expressed as a linear combination of the products of distinct ones. As- 
sume the contrary and consider monomials 1]^ . . . of degree iV — d — 1 which 
cannot be expressed as a linear combination of products with all distinct forms. 
Among those monomials we take a monomial m = Z^V . . having the maximal 
number of distinct forms in the product. Since m is not a product of all distinct 
forms, it should contain a form U in some power ki > 2. Given that kt > 2 
and the degree of m is — d — 1, one can find d + 2 distinct forms Z^^ , . . . , Z^^^^ 
that do not appear in m. Since any d -f- 2 of our forms span the dual space of 
M'^+i^ we can express U as a linear combination of li-^, . . . Now rewrite 

m as (ai • Zij H h ad+2 ■ hd+2) ^1' ■ • • 'f'""^ ■■■I'n ^ where ai ■ + ■ ■ ■ + ad+2 • 

= ^i- Thus we get an expression of m as a linear combination of monomials 
ttj ■ U-li^ . . . Zf*~^ . . . Z^, where each such monomial has more distinct forms than 
m. Each of such monomials can be expressed as a linear combination of products 
of all distinct forms, since m was chosen as a monomial with the maximal possible 
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number of distinct forms, which cannot be expressed in such a way. This is a con- 
tradiction. Therefore m can also be expressed as a hnear combination of products 
of all distinct forms. □ 



Below we consider the inverse problem for a weakly non-degenerate S, using 
notation from (|3.ip and p.3p . Here we no longer have a natural basis of all simplices 
sharing a common vertex v^r. Because of that we need to consider all N points and 
include one more linear form Zjv into the corresponding matrix L. Slightly abusing 
our notation, we denote by L the same matrix as before, although it contains one 
more (last) column corresponding to vyy. Similarly to notation (|3.3p . for a given 
set J of (i+ 1 linear forms we denote by L( J) the determinant of the corresponding 
(d -I- 1) X (d -I- l)-minor of L. 

We introduce the extended (^J^) x (^^j^) -matrix Mats with columns consisting 
of the coefficients of the homogeneous polynomial li-^{u) . . . ^i„_^_i(u) with respect 

to the monomial basis in the variables (uq, mi, . . . , Ud)- By Lemma 1191 Mats ha-s 
full rank, since it determines a surjective linear map onto HPol{N — d — l,d + 1). 
Thus Mats has a maximal minor with a non- vanishing determinant. Formula (j4.ip 
holds for the determinant of any maximal minor of Mats- 

Lemma 20. Let & be any set of (^J^) columns of Mats- We label each column 
T G 6 by the corresponding subset of the linear forms of cardinality 

N — d — 1- Then the determinant of the maximal minor Mats{&} formed by the 
columns of & is given by: 



det 



Matsie)] = k{&) ■ Yl 



(4.1) 



where k{&) is a constant (possibly equal to zero) depending only on the combina- 
torial structure of the [N — d — l)-tuples in the set 6. 



Proof- Fix the set © as above. In what follows we treat both sides of (j4.1l) as 
complex- valued polynomials in N ■ {d+1) variables, these variables being the entries 
of matrix L. 

We first show that every determinant L(J) divides det AIats{&) - Indeed, let 
J — {ji, . . . , jd+i} be a set of {d + 1) forms which has a nonempty intersection 
with any (N — d — l)-tuple of forms in (3. Let z = (zi, . . . , zjv.((j+i)) be a zero 
of the polynomial L(J), which means that forms Ij-^, - - . comprised of the 

corresponding coordinates of z are linearly dependent. Therefore, there is a non- 



zero vector Uq G 



such that Zjj^(uo) = ■•• = ^^^^^(uo) = 0. Consider the 



row vector (ug) consisting of ^) homogeneous monomials of degree N — d — 1 
evaluated at Uq. We notice that (uq) is in the kernel of Mats{&), as the product 
of (ug) with each column vector T G © of Mats{&) is equal to YijeT hi^o)] and 
every set T G © contains at least one of the forms Ij-^, - - - , Ij^^-^ in such a product. 



Thus det 



Mats{&) 



also vanishes at such z. 



We recall a well-known fact (cf. e.g. [IQ, Theorem 61.1]) that 



L( J) = det 



Mo / 

Ml 



Zji+N 



^2 

Zj2+N 



Ud \Zji+dN Zj^+dN 



Zjd+i + N 



Zjd+i+dN I 
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is an irreducible complex- valued polynomial in variables Zj-^ , Zj^ . . . , Zj^^^+jvci- Now 
if every zero of an irreducible polynomial . . . , zjv((;+i)) annihilates another 
polynomial ^(zi, . . . , 2jv((i+i)) then p divides q. We conclude that L(J) divides 

det Mats (6) ■ 

Using the fact that each L(J) is an irreducible polynomial and all L(J)'s are 
pairwise distinct (i.e. have distinct sets of projective zeros) we conclude that the 

product of L(J)'s in the right-hand side of (j4.ip divides det Mats{G) ■ 

Finally, the product of L(J)'s has the degree 



(d+1) 





N 










d+1 





VT e © TDJ^ 



(4.2) 



We observe that for each T G 6 the complementary set of d -I- 1 forms cannot 
be taken as a feasible J. We notice further that these complements are the only 
exceptions for the choice of J. Therefore as a feasible J we can pick any of (J^-^) 
{d+ l)-tuples except those (^J^) complements of a T G 6. Therefore, (j4.2p equals 



(d+1) 



N 
d+1 



N -1 
d 



The latter expression coincides with the degree of the polynomial det Mats{&) 

(assuming that it is not a zero), as AIats{6) has (^J^) columns and each entry is 
a homogeneous polynomial of degree N — d — 1 . 

Hence Mats{&) coincides with the product of L(J)'s up to a constant factor 
which might vanish. This constant does not depend on the entries of matrix L and 
hence it is completely determined by the set 6, regardless of the location of points 
of S in W^. □ 

Lemma [TO] allows us to solve the inverse moment problem for a given weakly non- 
degenerate S = {vi, . . . , vjm} in a certain linear space DJl{S) 3 DJl^{S) of measures 
supported on conv(S'). Namely, 9Jl(S') is spanned by measures fi G 971(5') whose 
normalized moment generating functions ^)i(u) belong to D\ai{S), i.e. -F),(u) = 
P(u)/<i>5(u), where $s(u) = IljLi^jlu) and P{u) is a polynomial of degree at 
most N — d — 1. Indeed, by Lemma [TOl anv -R(u) G fHot(S') can be represented in 
the form 



il<i2<---<id+l<N 
If V, , , V 



III ■ ■ 



ld + 1 



■ld+1 ■ 



'ld+1 



span 



(4.3) 



then the term 



with some real constants Ki-^^i^, 

— can be interpreted as the normalized moment generating function of an 



■ h 



appropriately scaled standard measure of the d-dimensional simplex spanned by 
these vertices. 

If Vij^ , , . . . , Vj^^j only span a hyperplane H in M"^ then (|4.3p corresponds 
to a singular (w.r.t. to the Lebesgue measure on M'*) measure fig supported on 
6 — conv(vij^, . . . , Vi^^j^). One way to define it as the weak limit of a sequence of 
(absolutely continuous with respect to the Lebesgue measure on W^) measures — the 
appropriately scaled standard measures /z^j of family of d-dimensional simlices St 
which degenerate into S when t = 0. There is no loss in generality in assuming 
Ki-ii2...id+i = Ij to deal with probability measures. 

Proposition 21. Let W = {wi, . . . , Wd, w^+i} be a (d + l)-tuple of points in M.'^ 
such that W spans a hyperplane H C M"*. Denote by = 1 ^ (wi, u), . . . , ^w^+i = 
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1 — (wtj+i, u) the associated linear forms. There exists a unique measure /iw sup- 
ported on S = conv(W) with the normahzed moment generating function (u) 
given by 

P^M = , , \ ■ (4.4) 

Proof. Without loss of generality assume that W = {wi, . . . , w^, Wd+i} is or- 
dered in such a way that {wi, . . . , w^} span H. Then, St is defined as St = 
conv(5, Wi^^j with z a unit normal to H, and fig^ as the uniform density proba- 
bility measure supported on St- Then lim fj,s^ = /i^, where lim is understood in sense 

of weak convergence of distributions (measures), i.e. that lim / /d/i^j = J fdfis for 

any bounded, continuous real function on R"*, cf. e.g. [9\. Then, this measure has 
compact support, and thus is determined by its moments, cf. e.g. [245 Proposi- 
tion 3.2]. □ 

Remark 12. One can prove that the integration of a smooth compactly supported 
function (p with respect to the limiting measure /iw is given by the integration of cf) 
over S with a continuous piecewise linear weight function uniquely determined by 
S. Similar limits appear frequently in the theory of splines. Since we only need the 
existence of /iw we do not pursue this topic here. 

Our solution of the inverse moment problem for the linear space 971(5') closely 
follows the pattern presented in Example [5] In other words, given a weakly non- 
degenerate S and the set of moments up to order N — d — 1 we 

(i) produce the rational function i?(u) S 9^ot(S') with Taylor coefficients coincid- 
ing with the normalized moments; 

(ii) represent i?(u) in the form (j4.3p : 

(iii) for each term as in (14. 4p determine the underlying measure supported on the 
(probably degenerate) convex hull of the vertices v^^ , v^^ , . . . , v^^^^ . 

We can now prove our central result claiming that D)l^{S) — Dyi{S) for a weakly 
non-degenerate S. 

Proof of Theorem [SI Theorem [5] is already settled in Corollary [TO] for the case of 
strongly non-degenerate S. It remains to consider the case of weakly non-degenerate 
S. The denominator of the moment generating function F-p{u) for an arbitrary 
generalized polytope V with the vertex set S is of the form H^iLih by Proposition [51 
and its numerator belongs to Pol{N — d — l,d). As S* is weakly non-degenerate, 
F-p(u) can be written as a linear combination of the fractions as in (|4.3p . where 
(ii, 12, . . . j(i+i) runs over the set of (d + l)-tuples of indices. If a (d -I- l)-tuple 
. . . is spanning then ^ ^ is the moment generating function of 

the measure supported on the simplex A, determined by its denominator, with the 
uniform density K/dWo\{A). By Proposition [211 if a (c? -f l)-tuple li-^^Ji^, . . .li^^^ 
is not spanning then — ; — —, is the moment generating function of a singular 

measure supported on a degenerate simplex. As P is a generalized polytope, its 
standard measure has no singular components. Therefore, no degenerate simplices 
can appear in its decomposition. □ 

Remark 13. The latter proof demonstrates that if one starts from the set of 
moments of the standard measure /x of a polytope with the vertex set S then we 
never obtain degenerate simplices while solving the inverse moment problem. This 
is why 971^(5') — D}l{S). However, an explicit description of 5"^ (5) for a general 
weakly non-degenerate S is missing at present. For concrete Examples [3] and |3] we 
give these descriptions below. 
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Our final result computes dim 971^(5') and describes a procedure to construct a 
basis for m'^iS). 

Proposition 22. Let S — {vi, . . . , v^} C Mf^ be an arbitrary weakly non-degene- 
rate spanning set. Then 

(i) dimOT'^(S') = {^-^) ~ ttdeg where fltjeg is the number of degenerate simplices, 
i.e. the number of non-spanning (d -I- l)-tuples of points of S. 

(ii) If 5 is a degenerate d-dimensional simplex with vertices in S \ {vi} then 
there is exactly one linear dependence among the standard measures of all 
d-dimensional simplices on and d vertices of S. 

(iii) The standard measure of any d-dimensional simplex on is contained in at 
most one dependence as in 

(iv) For any vertex one can construct a (in general, non-unique) basis Bi of 
9H^(S') consisting of standard measures of d-dimensional simplices on v^, as 
follows. 

(a) Start from the set Bi of the d-dimensional simplices on v^. 

(b) For each degenerate simplex S not containing v^, remove from Bi the 
standard measure of an arbitrary simplex on from the corresponding 
to 6 hnear dependence, cf. 

Thus we obtain {^^^) — ^deg standard measures of d-dimensional simplices, 
forming a basis of DJl^ (S) . 

Proof. To prove (0), notice that dimOT(S') = dim 9^01(5") = Cd^)- ^^11' 
971(5') = dyi{S) dJldeg{S), where dyideg{S) is the hnear span of the measures /i^^' 
with S running over the set of all degenerate simplices spanned by {d+ l)-tuples of 
dependent vertices in S, cf. Proposition Observe that these measures /i^^'' are 
linearly independent, as each degenerate simplex defines a singular measure sup- 
ported in a proper hyperplane, and these hyperplanes differ for different degenerate 
simplices. We are done with (P. 

Let Eq be a dependent d+ 1-subset of S, and S — conv(I]o) be as in ([n]). Then S 
spans a hyperplane Hq. As each d-dimensional simplex on do := and d vertices 
from So is uniquely defined by the latter, it suffices to analyze dependencies between 
the standard measures of d — 1-simplices with vertices in Eo- 

We can view Eo as a weakly non-degenerate subset in M"*^^ = Hq. By we 
have dim 971^ (So) = {^'^i) — tJrfeg(So). If So is strongly non-degenerate as a subset 
of = i.e. ^degi^o) — 0, then dim97t'^(Eo) = d, i.e. there is exactly one 

linear dependence between the standard measures of d — 1-simplices with vertices 
in So, and we are done. Otherwise, So = {(Ti}USi, with Si spanning a hyperplane 
Hi in Hq. Moreover, this can only happen if d > 3. Now, we can repeat the whole 
argument with ai in place of ctqi in place of Sq, and Hi in place of Hq. Again, we 
either have Si strongly degenerate, and we are done, or we repeat this argument, 
etc., until we hit a strongly non-degenerate S^, which is bound to happen, as the 
dimension goes down each iteration. This completes the proof of ([n]) . 

Then, (pli)) stems from the fact that the vertices of d-dimensional simplex on 
Vi distinct from V; span a hyperplane, and the only possibility for a degenerate 
simplex S as in ([u]) is to lie in this hyperplane. 

Finally, to prove (|iv]), observe that the set B'^ of the standard measures of 
d-dimensional simplices containing a given vertex always spans 97t^(S'), see 
Lemma [TSl Now for each degenerate d-simplex S we prune B'^ by removing the 
standard measure of a simplex in the linear dependence corresponding to S. In 
view of ([u]) and this process is well-defined and unambiguous. In the end we 
obtain (^J^) — tJdeg standard measures of d-dimensional simplices. In view of Q 
they form a basis of 97t'^ (S) , as claimed. □ 
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Remark 14. The above discussions show that the columns of Mats corresponding 
to degenerate simphces must necessarily be included in any non- vanishing maximal 
minor MatsiG). 

We conclude our discussion of the weakly non-degenerate case with two examples. 



V4 = (0, 2) V3 = (2, 2) V4 = (0, 2) 




V5 = (0,0) V2 = (2,0) vi=V5 = (0,0) V2 = (2,0) 

Figure 3. Vertices for Examples [3] and HI 



Example 3. Let S = {vi, V2, V3, V4, V5} where Vi = (1,1), V2 = (2, 0),V3 = 
(2, 2), V4 = (0, 2), V5 ^ (0, 0). Then = 1 - ui - M2, ^2 = 1 - 2wi, ^3 = 1 - 2mi - 
2u2,li = 1 — 2^2,^5 = 1. Calculating the products IJj, i < j and taking their 
coefficients in the standard monomial basis of Pol{2,2) we obtain the following 
6 X 10-matrix Mats 



Mats = 





hh 


^1^3 


hh 


hh 


^2/3 


hh 


hh 


hh 


^3^5 


hh 




1 


( 1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


\ 




-3 


-3 


-1 


-1 


-4 


-2 


-2 


-2 


-2 







U2 


-1 


-3 


-3 


-1 


-2 


-2 





-4 


-2 


-2 




ul 


2 


2 








4 



















UIU2 


2 


4 


2 





4 


4 





4 










ul 


I 


2 


2 














4 








/ 



Its rank equals 6 and one of non- vanishing maximal minors consists of the columns 
with numbers 6 = {5, 6, 7, 8, 9, 10}. (Recall that any non- vanishing maximal minor 
must include columns 6 and 9 corresponding to degenerate triples (vi,V3,V5) and 
(vi,V2,V4) resp.) The corresponding submatrix Mats{&) equals 



Mats{G) = 





hh 


hh 


hh 


hh 


hh 


hh 




1 


( 1 


1 


1 


1 


1 


1 


\ 




-4 


-2 


-2 


-2 


-2 







U2 


-2 


-2 





-4 


-2 


-2 




ul 


4 



















U1U2 


4 


4 





4 










ul 


I 








4 








/ 
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U1U2 u\ 
-4 4 \ 
-2 2 
-2 
-1 1 
-1 
-1 I 

Thus, given an arbitrary rational function R{ui,U2) = P(ui, W2)/4's(mi, U2) where 
P{u\,U2) = cioo + 011,0^1 + cio,iU2 + 012,0^1 + (iiiWi'U2 + O'02U2 is any polynomial of 
degree at most 2 and $s(ui,U2) = hhhUh, we obtain 

Wl45 = — aio — 020 — 0-11 

wi35 = -5 (an - 002) 

W134 = 5(001 - On) 

W125 = -3(000 - oio - On + 002) 

■u;i24 = 2(010 - On) 

,Wl23 = -3(010 - Ool - 020 + On). 

Triangles A135 and A124 are degenerate which implies that if the original measure 

we arc recovering is polygonal then W135 = W124 = 0. Therefore the linear space of 
numerators P{ui,U2) for the space in this example is given by the relation 

ooi = On = oo2- 

Our last example is more degenerate! than the previous one, although still weakly 
non-degenerate. In fact, in this example 5 is a multiset since Vi = V5. It shows 
that our technique can be generalized to a certain class of multisets as well. 

Example 4. Let S = {vi, V2, V3, V4, V5} where vi = V5 = (0,0), V2 = (2,0), V3 = 
(1, 1), V4 = (0, 2). Then = = 1, ^2 = 1 - 2ui, ^3 = 1 - ui - U2, ^4 = 1 - 2u2. 
Calculating all products lilj, i < j and taking their coefficients in the standard 

monomial basis of Po/(2, 2), we obtain the following 6 x 10-matrix Mats'- 

^3^4 ^3^5 ^4^5 

1 1 1 \ 
-1 -1 
-3 -1 -2 

0" 

2 
2 / 

Its rank equals 6 and a non- vanishing maximal minor consists of the columns with 
numbers & = {1,3,4,5,6,8}. The corresponding submatrix Mats{&) equals 

^2^3 ^2^4 ^3^4 

1 1 1 \ 
-3 -2 -1 
-1 -2 -3 

2 
2 4 2 
2 / 



Further, 



AMats (S) 
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Ml 


U2 


ul 


/2/3 


/o 
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Z2Z4 














hh 








2 





I3I4 


1 


-1 








hh 





1 








hh 




-1 


1 


1 



Mats = 





hh 


^1^3 


hh 


hh 


^2^3 


hh 


hh 
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( 1 
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-1 








-3 


-2 





U2 
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-2 





-1 


-2 


-2 
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U1U2 
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Further, 
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^? 


UiU2 






hh 
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hh 
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-2 




hh 
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hh 
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hh 











-1 


1 


-1 




hh 


Vo 
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/ 



4Matg (6) 



Thus, given an arbitrary rational function R(ui, U2) = P{ui, U2) /$5(ui, M2) where 
P{ui,U2) = ooo + cLifiUi + flQ 1U2 + a2.ouf + aiiUiU2 + ao2u\ is any polynomial of 
degree at most 2 and $s(mi,M2) = hhhhh, we obtain 



u'345 = -3(2010 
^^235 = -i((2aoi 
W234 = 4aoo 
W145 = 2a2o 

1^125 = 2ao2- 



2a2o + an) 
' an + 2ao2) 
2aoi + 2aio + 020 + an + ao2 

an — 020 



Notice that triangles A125, A135, A145, A234 are degenerate. If we know that the 
original measure we are recovering is polygonal then one should get wi2b = ""^135 = 
tfi45 = W234 = 0. Therefore, the linear space of numerators P{ui,U2) for the space 
S^{S) in this example is given by the system of equations: 




3aio — 0. 



5. Remarks and open problems 

Remark 15. A weaker form of Corollary U] (i.e. the rationality of F^{u), but 
without the claim on the particular shape of the denominator) can be derived 
directly from (|1.3p by using Stokes formula, along the lines of [6, Lemma 1]. 

Problem 2. Find an appropriate version of Theorem [51 applicable to non-simple 
and/or non-convex polytopes. 

Remark 16. Choose an arbitrary basis {A^} of Dyi^{S) consisting of the standard 
measures of simplices. The set {Aj} spans an integer lattice in Dyi^{S). (One can 
easily see that this lattice is invariantly defined independently of the choice of a 
basis of standard measures of simplices.) Denote by DJl^{S) the space 971^(5) with 
the latter lattice. We can prove the following. 

Proposition 23. If Conjecture [7] holds then any generalized polytope V G 7^(5') 
with standard measure /ip in OT^(S') corresponds to a rational point in 9Jt^(S'). 

Proof. One can easily show that any P as above can be represented as the union of 
the closures of connected components of R"^ \ H{S), where H{S) is the hyperplane 
arrangement consisting of all hyperplanes spanned by d-tuples of points in S. (The 
converse is obviously not true.) Let S ^ S he the extended set of vertices obtained 
by adding to S all vertices of the hyperplane arrangement H{S). Since each con- 
nected component in R'^ \ H{S) is convex, it can be triangulated on S. Consider 
the space DJl^{S). Obviously, fi-p is an integer point in DJl^{S). Also, 9Jl^{S) is 
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contained in dJl^{S) as a sublattice. Thus, if fi-p belongs to DJl^{S) it is a rational 
point there. □ 

Remark 17. If either V itself or its complement in conv(S') can be triangulated 
by simplices with vertices in S then ii-p is an integer point in DJtz{S). 

Let us pose the following tantalizing question. 

Problem 3. Does there exist a generalized polytope V with the set of vertices 
S{V) and with non-integer coordinates in dJlz{S{V))7 

Problem 4. One can also define a rational convex cone ^os{S) C dJlz{S) by 
taking non- negative linear combinations of all fx-p, where V runs over the set of all 
generalized polytopes in 'P{S). 

Conjecture 24. The rational cone *Pos(S') is uniquely determined by the oriented 
matroid associated to S. 

We conclude this section with the following question. 

Problem 5. Is it possible to describe the extremal rays of ^os{S)l 

One can easily show that a simplex from 7^(5) spans an extremal ray of Vjios{S) if 
and only if it does not contain any points of S distinct from its vertices. Problem[5] 
is apparently closely related to the problem of classification of combinatorial types 
of point arrangements, see e.g. |14) and references therein. 
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